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In this work, we show that an in-plane magnetic field can drive a fully gapped p±ip topological 
superconductor into a gapless phase which supports symmetry protected Majorana edge states 
(MESs). Specifically, an in-plane magnetic field can close the bulk gap and create zero energy 
Majorana flat bands (MFBs) in the excitation spectrum. We show that the MFBs in the gapless 
regime are protected by a chiral symmetry and are associated with MESs. Novel uni-directional 
MESs which propagate in the same direction on opposite edges appear when the chiral symmetry 
is broken. The MFBs and the uni-directional MESs induce nearly quantized zero bias conductance 
in tunneling experiments which are robust in the presence of a gapless bulk and disorder. 



Introduction — A topological superconductor (TS) has 
a bulk superconducting gap and topologically protected 
gapless boundary states [IHS] . TSs are under intense the- 
oretical and experimental studies due to the possibility of 
realizing Majorana fermions in these systems, which act 
as their own antiparticles and obey non-Abelian statis- 
tics [6-8 . Majorana fermions in TSs are topologically 
protected, in the sense that the Majorana fermions can- 
not be removed by perturbations unless the bulk energy 
gap is closed or certain symmetries are broken. 

Remarkably, recent development shows that topologi- 
cally protected Majorana fermions exist in systems where 
the bulk is gapless [9tiT9] . For example, Majorana 
edge states (MESs) with flat dispersion can be found 
in 2D nodal d xy + p-wave superconductors which re- 



a) 



spect time-reversal symmetry [TQ1TL2] . It is also shown 
that zero energy Majorana flat bands (MFBs) can ap- 
pear on the surface of 3D time-reversal invariant non- 
centrosymmetric superconductors which have topologi- 
cally stable line nodes in the bulk [T3HT5] . Proposals on 
gapless TSs which break time-reversal symmetry have 
been made recently [T6HT8] . 

All the cases mentioned above are intrinsic gapless su- 
perconductors. In this work, we show that an in-plane 
magnetic field can drive a fully gapped p ± ip-wave TS 
into a gapless regime which supports symmetry protected 
Majorana fermions. An in-plane magnetic field may first 
close the bulk gap. Further increasing the strength of 
the magnetic field creates zero energy MFBs in the exci- 
tation spectrum when the bulk is gapless. Importantly, 
we show that the MFBs are protected by a chiral sym- 
metry and are robust against disorder. The evolution of 
the excitation spectrum of a p ± ip-wave superconductor 
as a function of the in-plane magnetic field strength is 
shown in Fig.l and Fig. 2. 

S-wave pairing and Rashba spin-orbit coupling terms 
break the chiral symmetry which protects the MFBs and 
lift the Majorana fermions from zero energy. In this case, 
uni-directional MESs, which are distinct from the usual 
helical or chiral MESs in that the modes on opposite 
edges move in the same direction, may appear. Finally, 
we show that the MESs survive in the presence of disor- 
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FIG. 1: a) A schematic picture of a p ± ip-wave supercon- 
ductor subject to an in-plane magnetic field V y . A tunnel 
junction and a normal lead N is attached to the superconduc- 
tor, b) The energy spectrum of a p zb ip superconductor in 
the topologically non-trivial regime. Periodic boundary con- 
ditions in the x-direction and open boundary conditions in 
the y-direction are assumed. The parameters are t = 12 A p , 
ji — 3A P — 2t, A s — 0, cxr — and V y — 0. c) Same parame- 



ters as b), except V y 
in this regime. 



0.7A P . The bulk energy gap is closed 



der and induce nearly quantized zero bias conductance 
in tunneling experiments. 

Majorana Flat Bands — We start with a BdG Hamil- 
tonian which describes ap± ip-waxe superconductor 



H p (k) 
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Here £(k) = [— t(cosk x + cosk y ) — /j]ao — 
sin k y a x -f sin k x a y ] is the sum of the kinetic energy 
and the Rashba spin-orbit coupling, V describes the Zee- 
man coupling of the electrons with an external magnetic 
field, A(k) = (A s + d(k) • a)(ia y ) is the superconduct- 
ing gap function. We first assume that the spin-singlet 
pairing amplitude A s and the Rashba spin-orbit coupling 
ctR are zero. The spin-triplet pairing vector is chosen as 
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FIG. 2: The evolution of the energy spectrum of a p±ip-w&ve 
superconductor as V y increases. For a) to c) the parameters 



are the same as Fig. lb except the values of V y . 



Vy A v 



b) V y = 3A p . c) V y = 3.5A P . d) S-wave pairing and Rashba 
terms with A s = 0.3A P and cxr = 0.2A P are added to a). 



d(k) = A p (— sin k y , sin k x , 0) such that the Hamiltonian 
descirbes a two dimensional p ± ip-waxe superconductor 
where A p is a constant. When V = 0, the Hamilto- 
nian respects both time-reversal symmetry T = UtK 
with /7~ 1 i^*(k)/7 T = H p (—k) and particle-hole symme- 
try P = U P K with U^H^Up = -H p (-k). Here, K 
is the complex conjugate operator, Ut = &o <£> io~ y and 
Up = cf x (g) cr such that T 2 = —1 and P 2 = 1. 

According to symmetry classification, the above 
Hamiltonian in the absence of an external magnetic field 
belongs to Dili class which can be topologically non- 
trivial pp. In the topologically non-trivial regime, the 
p ± ip superconductor possesses gapless counter propa- 
gating helical MESs. The energy spectrum in the topo- 
logically non-trivial regime is shown in Fig. lb. In the rest 
of this section, we show that the p ± ip superconductor 
responds to an in-plane magnetic field in an anomalous 
way as described in the Introduction. 

To be specific, we suppose a magnetic field is ap- 
plied in the y-direction such that V = (0, V y ,0). In the 
presence of a magnetic field, the time-reversal symme- 
try T = UtK is broken. However, one can show that 
the Hamiltonian satisfies a time-reversal like symmetry 
T ld = U T idK such that T^H{k x , k y )T ld = H(k x: -k y ), 
where Urid = cr z <&cr z . Moreover, the Hamiltonian satis- 
fies a particle-hole like symmetry P\ d = U p\ d K such that 
P~^H(k x ,k y )Pi d = -H{k x ,-k y ) with U P i d = cr y ®a y . 
Due to the fact that the symmetry operators operate on 
k y only and k x is unchanged, one may regard k x as a 
tuning parameter and the Hamiltonian can be written 
as Hk x (k y ). Since Hk x (k y ) satifies the symmetries T\ d 
and Pi d with T 2 d = P 2 d = 1, H kx (k y ) is a BDI class 
Hamiltonian which can be classified by an integer pp. 

To classify the Hamiltonian Hk x (k y ) with k x as a tun- 



ing parameter, we note that as a result of the T\ d and 
Pid symmetries, Hk x (k y ) satisfies the chiral symmetry 
Sid = T ld P ld with 



$i d H(k x , ky)S\ d — H(k x ,ky). 



(2) 



In this case, Hk x (k y ) can be off-diagonalized in the basis 
which diagonalizes S\ d such that 



HkAky) -[Al(k v ) 



(3) 



Note that Ak x (k y ) is real at k y = 0, ±7r, we can define 
the quantity 

z(k) = = Bet[A kx (k)]/\Bet[A kx (k)}l (4) 

such that 6{k) = nn at k = 0, ±7r with integer n. The 
winding number of 6(k) can be used as the topological 
invariant which characterizes the Hamiltonian Hk x (k y ). 
The winding number Nbdi can be written as [20] 



Nhdt = — 



i f k y^ dz(k y ) 



k Jk v =o z (ky) 



(5) 



Using A kx (k y ) obtained from H kx (k y ), it can be shown 
that \Nbdi\ = 1 when 



M(k x , k y — 0)M(k x , k y = tt) < 0, where 
M(k x ,k y ) = [/i + t(cos k x + cos ky)} 2 1 A ' 2 



A 2 sin^ k x 



(6) 

assuming that V y and A p are non-zero. In the range of 
k x where Nbdi = 1, the Hamiltonian H kx (k y ) is topo- 
logically non-trivial. For a p ± ip superconductor with 
periodic boundary conditions in the ^-direction and open 
boundary conditions in the ^/-direction, there are zero en- 
ergy Major ana modes localized on the edges of the system 
when Eqj6] is satisfied. Therefore, MFBs appear in the 
corresponding parameter regime. 

The evolution of the energy spectrum of a p ± ip su- 
perconductor as a result of an increasing in-plane mag- 
netic field is shown in Fig.l and Fig. 2. First, an in-plane 
magnetic field reduces the bulk gap as shown in Fig.lc. 
Second, after the bulk gap is closed, MFBs appear for a 
finite range of k x where \Nbdi\ = 1 as shown in Fig. 2a. 
Third, by further increasing the magnetic field, the bulk 
gap at k x = is closed (Fig. 2b). Fourth, by increasing 
the magnetic field even further, the energy crossing at 
k x = disappears and only a MFB remains (Fig. 2c). It 
is important to note that the MFBs appear when the 
bulk is gapless. The bulk energy spectrum of a p ± ip- 
wave superconductor corresponding to Fig. 2a is shown in 
Fig. 3a. It is evident that there are nodal points in the 
bulk spectrum when MFBs appear. The nodal points in 
Fig. 2a are the projection of the bulk nodal points on the 
/c^-axis in Fig. 3a, similar to the cases in intrinsic gap- 
less TSs [H US [19]. Both the nodal points in the bulk 
spectrum as well as the MFBs are protected by the topo- 
logical invariant Nbdi- In other words, the MFBs and 
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FIG. 3: a) The bulk energy spectrum of a p ± ip supercon- 
ductor. The parameters are the same as the ones in Fig. 2a 
but periodic boundary conditions in both the x and y di- 
rections are imposed, b) The energy spectrum of a pure 
p =b ip-wave superconductor with dimensions L x = 100a and 
L y = 300a, where a is the lattice spacing. Only the thirty 
lowest energy eigenvalues are shown. Periodic boundary con- 
ditions in the x-direction and open boundary conditions in 
the ^/-direction are assumed. On-site Gaussian disorder with 
variance w 2 = (1.5A P ) 2 is present. It is evident that as Vy 
increases, states collapse to zero energy and stay there, in- 
creasing the number of zero energy modes. This indicates the 
widening of the MFBs as V y increases. Importantly, the zero 
energy Majorana modes are not lifted by disorder. 

the nodal points in the bulk appear whenever Nbdi is 
non-trivial for some range of k x . 

Uni- directional Majorana Edge States — It is shown 
above that MFBs appear when Nbdi = 1 for a finite 
range of k x and the MFBs are protected by the chiral 
symmetry Sid in Eqj2] However, s-wave pairing and 
Rashba terms, which can arise in non-centrosymmetric 
superconductors [2T] . break the chiral symmetry and 
lift the zero energy modes to finite energy as shown 
in Fig. 2d. In the case of adding s-wave and Rashba 
terms to Fig. 2c, the MFB acquires a finite slope and uni- 
directional MESs appear at the sample edge as shown in 
Fig. 4a. A schematic picture of the uni-directional MESs 
is shown in the insert. 

Another interesting finding is that the uni-directional 
MESs can appear in the absence of p ± ip-wave pairing. 
The energy spectrum of an s-wave superconductor with 
Rashba terms and finite V y is shown in Fig. 4b. It is ev- 
ident that uni-directional MESs appear in this case. To 
understand the origin of the MESs in the gapless phase, 
we note that the Hamiltonian H p (k x ,k y ) at k x — sat- 
isfies the chiral symmetry S = PT\& and is classified by 
the topological invariant Nbdi when A p = 0. If both 
the A s and qlr are non-zero, the Hamiltonian at k x = 
has winding number \Nbdi\ = 1 when 

M 8 (0, 0)M 3 (0, tt) < 0, where , . 

M s (k x , k y ) = [/x + t(cosfc a . + cos k y )] 2 + A 2 s -V 2 . t ' 

Non-trivial Nbdi indicates the appearance of the uni- 
directional MESs. It is important to note that sys- 
tems with pure s-wave pairing and Rashba terms can be 




FIG. 4: a) The parameters are the same as those in Fig. 2d 
except the values of V y . The MFB acquires a finite slope when 
A s and otR are finite at V y = 4A P . Uni-directional MESs ap- 
pear in this regime. A schematic picture of the uni-directional 
MESs is shown in the insert, b) The energy spectrum in the 
presence of spin-singlet pairing and Rashba and no p-wave 
pairing terms. Here, A p = 0, t = 12 A s , \i — 3A S — 2t and 
a R = 2A s . 



realized by inducing s-wave superconductivity in semi- 
conductors as demonstrated in recent experiments [22]- 
[24] . This opens a way for realizing the novel uni- 
directional MESs. 

Experimental Detection and Effects of Disorder — 
It has been shown in previous works that Majorana 
fermions induce resonant Andreev reflection at the junc- 
tion between a normal lead and a fully gapped TS 
[25] GS]- However, resonant Andreev reflection may not 
happen when the bulk is gapless due to the non- vanishing 
direct tunneling amplitudes from the normal lead to the 
gapless superconductor. Finite direct tunneling ampli- 
tudes make the reflection matrix non-unitary and render 
the arguments leading to resonant Andreev reflection not 
applicable [26]. In this section, we calculate the zero 
bias conductance (ZBC) of a junction between a normal 
lead and a TS as a function of the in-plane magnetic 
field strength. It is found that MFBs and uni-directional 
MESs induce nearly quantized ZBC even when the bulk 
is gapless and in the presence of disorder. 

A schematic picture of the experimental setup is shown 
in Fig. la. A normal lead is coupled to an edge of the TS 
to form a NS junction. Using the lattice Green's func- 
tion method [27H29] . we calculate the direct tunneling 
amplitude and the Andreev reflection amplitude of the 
NS junction. The results for the ZBC, ^ at zero volt- 
age bias V = 0, are depicted in Fig. 5. 

Fig. 5a shows the ZBC as a function of V y for a p ± ip- 
wave superconductor in the presence of on-site disorder. 
The ZBC depends on V y and the number of channels N c 
in the normal lead. The ZBC for a lead with N c = 4 (in- 
cluding spin degeneracy) is plotted in Fig. 5a. The direct 
tunneling contribution to the ZBC, T, is also shown. 

To understand the results, we note that time-reversal 
symmetry is preserved and the system is fully gapped at 
V y = 0, the 4e 2 /h quantization of ZBC is the property 
of a Dili class TS which has two Majorana zero modes 
on the edge [30, 31 . As V y increases, time-reversal sym- 
metry is broken and the ZBC is suppressed by disorder. 
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FIG. 5: The ZBC versus V y . a) ZBC from a lead to a pzLip- 
wave superconductor. The setup is depicted in Fig. la. The 
superconductor has dimensions L x = 100a and L y = 300a. 
Periodic boundary conditions in the x-direction is assumed. 
t = 12 A p , \i — 3A P — 2t in the superconductor and the lead. 
The barrier between the lead and the superconductor is sim- 
ulated by a reduced hopping amplitude t c = 0.3t. A semi- 
infinite lead with width 8a is used in the simulation. The 
number of conducting channels in the lead is N c — 4. The di- 
rect tunneling contribution to the ZBC is denoted as T. Gaus- 
sian on-site disorder with variance w 2 = (1.5A P ) 2 is present, 
b) S-wave pairing and Rashba terms with A s = 0.3A P and 
aR = 0.2A P are added to a). The red lines near V y = 4A P 
indicate the regime where the ZBC is quantized due to the 
presence of the uni-directional MESs. 

However, further increasing V y closes the bulk gap. When 
this happens, there is a large jump in the ZBC. This jump 
is due to the contribution from the Andreev reflection 
caused by the MFBs and the direct tunneling caused by 
the gapless bulk. This can be clearly seen from the V y 
dependence of T in Fig. 5a. It is interesting to note that 
the final ZBC is almost quantized at with N c = 4 

due to the Andreev reflection caused by the large number 
of independent Major ana fermions from the flat band. In 

9 2 

general, the ZBC caused by MFBs saturates at ^-N c . 

The nearly quantized ZBC at large V y in Fig. 5a sug- 
gest that the MFBs are robust against disorder. To con- 
firm this, the energy spectrum of the p±ip superconduc- 
tor with parameters corresponding to Fig. 5a is shown 
in Fig. 3b. It is evident from Fig. 3b that finite energy 
states collapse to zero energy and stay there, increasing 
the number of zero energy modes as V y increases. 

Indeed, the MFBs are protected by the chiral symme- 
try S\d m Eqj2] It can be shown that zero energy edge 
states, which are simultaneous eigenstates of the Hamil- 
tonian and Sid, have the same chirality for states local- 
ized on the same edge of the sample. The chirality of 



an eigenstate of Sid is defined as the eigenvalue of the 
state on Sid, which is always +1 or —1. The net chiral- 
ity number (NCN) of an edge of the sample, which is the 
sum of the chirality numbers of all the zero energy states 
localized on the edge, is always non-zero when MFB ap- 
pears. Moreover, it can be shown that the number of 
stable zero energy modes on an edge equals the NCN of 
the edge [12 j. Since on-site disorder does not break the 
chiral symmetry and cannot change the NCN, the num- 
ber of stable zero energy modes cannot be changed by 
disorder. 

It can be shown in the p±ip superconductor case that 
opposite edges have opposite NCNs. Therefore, the zero 
energy modes can be removed only if the chiral symme- 
try Sid is broken or when two edge states with opposite 
chiralities are coupled to each other. A similar index 
theorem for time-reversal invariant TSs is first shown in 
Ref . [T2] . However, in the time-reversal invariant case, 
the NCN on an edge is always zero and the flat bands 
are not robust against disorder. 

Fig. 5b shows the ZBC versus V y when s-wave pairing 
and Rashba terms are added to the Fig. 5a. At V y =0, the 

ZBC is quantized at 2^- as expected for a fully gapped 
Dili class TS [301131]. When V y is increased, the ZBC 
decreases due to disorder. Further increasing V y closes 
the bulk gap and there is a jump in the ZBC. Moreover, 
it is important to note that there are ZBC plateaus near 
V y = 4A p . As it is shown before, uni-directional MESs 
appear near this regime as shown in Fig. 4a. Even though 
the bulk is gapless, the number of bulk states with zero 
energy is very small for certain parameter regimes. In 
this case, the direct tunneling amplitude is small and 
one can ignore direct tunnelings. As a result, the ZBC 
is quantized due to Majorana fermion induced resonant 
Andreev reflection [25] . 

Conclusion — We show that an in-plane magnetic field 
can drive a p±ip-wdive superconductor to a gapless phase 
which supports MFBs. In the presence of s-wave pairing 
and Rashba terms, the MFBs acquire finite slopes and 
uni-directional MESs appear. These Majorana modes are 
symmetry protected and robust against disorder. They 
induce nearly quantized ZBC in tunneling experiments. 
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